Tunneling edges at strong disorder 
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Scattering between edge states that bound one-dimensional domains of opposite potential or flux is 
studied, in the presence of strong potential or flux disorder. A mobility edge is found as a function 
of disorder and energy, and we have characterized the extended regime. In the presence of flux 
and/or potential disorder, the localization length scales exponentially with the width of the barrier. 
We discuss implications for the random-flux problem. 

PACS numbers: 73.20.Jc,73.40.H 



The remarkable transport properties of a 2d electron 
gas in a strong perpendicular magnetic field are well- 
known [0: jumps in Hall resistance and peaks in lon- 
gitudinal resistance occur as the Fermi level crosses the 
single energy, at the center of the Landau band, where 
extended states are believed to lie. 

Experiments and numerical studies show that, as the 
Fermi level Ec + 6E approaches the center of the band, 
Ec, the correlation length C diverges: ( oc \SE\~'^'', where 

~ 7/3. Disorder plays a crucial role in all models for 
this transition. In a strong magnetic field B, eigenstates 
are confined to equal-energy contours of a disordered po- 
tential, slowly-varying on scales of order the magnetic 
length Ib = ific/eBf/'^. For a generic random, sym- 
metric potential V{r) with characteristic length |Q Aq, 
each of these equipotentials, for energies not too close to 
Ec = (T^), will traverse a region of finite size, exceeding 
the cluster diameter a(|(5i?|) in cross-section with expo- 
nential rarity. When 6E —> 0, long-range spatial correla- 
tions emerge in the equipotentials as the cluster diameter 
diverges at the continuum percolation transition Q . 

A model for this transition was proposed by Trug- 
man who argued that one might neglect tunneling 
among distinct equipotentials because of the small over- 
lap of their respective wave- functions. If this omission 
were valid, then the localization length would diverge as 
the spatial extent of the percolating eigenfunctions. The 
correlation exponent Vc for the classical percolation tran- 
sition is known to be 4/3. 

Mil'nikov and Sokolov [MS] attempted to calcu- 
late the correction to Trugman's classical picture arising 
from quantum tunneling between equipotentials. They 
argued that the dominant contribution of tunneling de- 
rived from the saddle-points of the potential along the 
largest eigenstates. Since these saddle-points occur be- 
tween longest equipotentials that traverse distances of or- 
der a{\6E\) = \6E\~'^^'^, and the correction they obtained 
from tunneling was inversely proportional to \SE\, MS 
suggested that 7/3 = 4/3+1. 

Their proposal, even were it correct, is not sup- 



ported by their argument. While saddle-points between 
longest equipotentials are rare, each equipotential is also 
a (bond) percolation hull, bounding a complex, fractal 
network of wave- functions with energy Ec- MS offered 
no justification for their neglect of the scattering between 
the hull and this internal network, which could lead to 
localization via back-scattering. 

The most widely accepted model for the Hall transi- 
tion, the Chalker network model , yields numerically a 
i/q close to 7/3. The Chalker model consists of a regular 
lattice of tunneling saddle-points, to be contrasted with 
the 'topologically disordered' network of equi-contours 
that underlies the classical percolation transition. This 
contrast with classical percolation has lead to the claim 
that quantum tunneling fundamentally alters the classi- 
cal transition, and that the relation Vq ^ 1 + Vc is merely 
coincidence. 

This issue assumes special importance for understand- 
ing eigenstates in a random-flux background (with van- 
ishing mean flux). If quasi-classical arguments applied, 
one would expect an extended state in the presence of 
random-flux that exhibited much the same properties 
as obtained for potential disorder; however, with appro- 
priate modification to incorporate the vanishing net chi- 
rality of the random-flux problem, the Chalker model 
displays no extended states 

Thus, the Chalker model appears to suggest that a pic- 
ture of of wave- functions, whose spatial extent is domi- 
nated by classical equipotentials, is naive. This state of 
affairs prompted us to reconsider tunneling among dis- 
tinct equipotentials. In particular, we would like to un- 
derstand when it may be possible to think of equipoten- 
tials as distinct classical units, and to calculate trans- 
port properties of a network of eigenstates by perturbing 
around their spatial configurations. 

Such an understanding might be relevant to a vari- 
ety of physical systems: not only to the QHE and the 
random-flux problem, but also to gauge-field theories 
of flux phases, where boundaries between distinct flux 
domains create a network of one-dimensional edges 
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Edges created by nonhomogeneous magnetic field in the 
absence of disorder or in the context of the quantum 
Hall effect with small disorder ||] have been studied else- 
where; here we study tunneling between equipotentials 
in the presence of strong potential and/or flux disorder. 

We construct a pair of one-dimensional wave-functions 
on a strip of length L and width M lattice constants a. 
Our lattice Hamiltonian is: 

H = Y,V^\i){i\+{J2u,\i){j\+h.c.} (1) 

where Vi denotes the on-site energy, and the hopping 
matrix elements Uj vanish except for nearest neighbors. 
Finite magnetic flux is included by allowing complex tij , 
where the local phase of Uj is determined by the flux en- 
closed within a square plaquette of the lattice. (Energies 
are in units of \t\ = 1, and we also take the lattice con- 
stant a = 1). We calculate localization lengths for this 
system numerically by computing the largest eigenvalue 
of a product of L transfer matrices, each of which adds a 
single row of M lattice sites at the end of the strip, and 
finite-size scaling j^]. These methods are now standard, 
although to our knowledge they have not been applied 
before to the geometries we discuss here. 

A pair of parallel edges, separated by roughly AI/2 
lattice sites (one edge centered between sites M/2 and 
AI/2 + 1, the other between sites M and 1) was con- 
structed in several ways {i denotes transverse coordinate 
on the strip): 

(la) Two parallel stripes of flux Bq with equal magni- 
tude and opposite sign: Bi — —Bq9{M /2 — i) + BQ9{i — 
M/2), random potential V is chosen on the interval 
[-W/2,W/2]. 

(16) Two parallel stripes of uniform potential Vq and 
width M/2, with equal magnitude and opposite sign in 
uniform flux Bq. Superposed on the potential stripes is 
random potential, chosen from [—W/2, W/2]. 

(2a) Two parallel stripes of flux, with flux through 
plaquette i chosen randomly from [0, Bq] for i < M/2, 
and from [-Bq, 0] for i > M/2. 

(26) Two parallel stripes of uniform potential Vq and 
width equal in magnitude and opposite in sign. Super- 
posed on the potential stripes was random flux with finite 
mean Bq, chosen from the interval [Bq — 5B,Bq + 5B]; 
\5B\ < \Bo\ 

All random quantities were chosen uniformly on the 
specified intervals, and independently on each lattice 
site. Periodic boundary conditions were imposed, except 
where explicitly noted below. We have concentrated in 
the regime of strong disorder, with W on order of the 
inter-Landau subband spacing. 

Flux and potential configurations differing locally from 
those listed above over distances of a few lattice constants 
did not affect the scalings we discuss below. We also al- 
lowed our edges to meander randomly in the y direction. 



and their transverse separation to fluctuate. Provided 
that the y coordinate of an edge was allowed to vary by 
at most a few lattice constants for each lattice site tra- 
versed in the longitudinal direction, our scalings similarly 
remained unchanged from those we found for configura- 
tions. We emphasize that our results are insensitive to 
details. The features of our edges that are important for 
the following discussion appear to be their low curvature, 
and their local (mean) flux/potential gradient, which re- 
mains constant as we increase the width of the barriers. 

We studied the dependence of the longitudinal local- 
ization length upon S, the width of the stripes in units 
of the lattice constant. Although we have for simplicity 
described geometries with two stripes only [S = M/2), 
our results are unchanged for multiple pairs of stripes. 

FIG. 1. Localization length versus energy for W — 5.0 (a) 
and W = 5.5 (b), and different values of strip width; M = 8 
(open circles), 16 (full circles), 32(squares), 64 (diamonds). 
The localization length Am has a relative error of 2% 

Figure 1(a) displays, for geometry (la). Am as a func- 
tion of E and stripe width. The intersections of the 
curves indicate mobility edges: in the localized regime, 
\m c)c Aoo, a constant. For Bq — 7r/4, this mobility 
edge disappears as indicated by figure 1 (b) when the po- 
tential disorder W exceeds roughly ~ 5.5, yielding 
localization at all energies. Figure 1(a) suggests that in 
the extended regime, Xm on M" with a ~ 2; however, 
we found that for larger M , Am increases exponentially, 
scaling as Am csc M"exp(M/A), where A represents a 
length scale determined by properties of the bulk. This 
exponential increase is apparent also for smaller values 
of disorder {W <C Wcr)- Scaled curves displaying this 
behavior are shown in figure 2. 

The critical properties of this transition will be dis- 
cussed elsewhere |l0| ; our present concern is the scaling 
with stripe width in the extended regime only. 

FIG. 2. Scaled localization length in the extended regime 
for W = 2.0, and different energies: E — 1.45 (full squares), 
3. (open circles), 1.475 (diamonds), 1.5 (triangles), 2.8 (open 
squares), 2.6 (crosses). Curves were collapsed by successive 
displacement along the horizontal axis; at each step the inter- 
polated (not extrapolated) displacements were minimized by 
least-squares. 

The scaling we find for M < A differs from that ob- 
tained for extended bulk states in potential disorder and 
uniform fiux with periodic boundaries, where Am oc M 
is expected at a band center, and is asymptotically con- 
stant elsewhere. (This latter geometry, in contrast to the 
geometries studied here, is isotropic). A Xm/M increas- 
ing with M above the band center has been observed 
previously for uniform fiux, Dirichlet boundary condi- 
tions, and potential disorder; however, those studies at- 
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tempted no quantitative analysis of this phenomenon Q . 
We have found numerically that this geometry also yields 
M"exp(M/A) scaling, and that this behavior is consis- 
tent with the earlier published data ^ . 

We obtained this exponential scaling for all edge con- 
figurations we studied, including (la, 6) and (2a, 5). Vq 
was set to 1.0 or 2.0 in cases (lb), {2b); the values of 
Bo and Vq were unimportant provided that the Landau 
subbands were not separated by much less than the char- 
acteristic scale of variation of W. 

When boundary conditions of the form la-26 or Dirich- 
let are imposed, chiral states traverse the boundaries. 
In the absence of edge-edge scattering, extended states 
propagate along these edges; the only mechanism for lo- 
calization of the edge wave-functions is scattering with 
states of opposite chirality at the other edge. We now 
argue that this inter-edge hopping gives rise to the scal- 
ing we have observed. 

At the quasi-classical level, an edge state lives on a 
line of constant potential. For strong disorder, an equipo- 
tential can wander into the bulk and meet the opposite 
edge, yielding backscattering. On an infinite plane, these 
equipotentials (percolation hulls) circumscribe islands of 
finite extent, the typical size given by SE"^ , SE the energy 
of the state with respect to the center of a magnetic sub- 
band. At the center of a subband, there exists an infinite 
equipotential that percolates through the sample. In this 
picture, the edge-to-edge scattering will depend on the 
probability of finding islands of size larger than S. The 
exponential rarity of these islands, sufficiently far from 
the percolation threshold, accounts for the exponential 
increase in localization length. If the typical island size £ 
is i ^ S, then we expect Ae to display algebraic behavior. 

When the disorder is small enough, ^ A, A the 
gap between magnetic subbands, edge states of energy 
Eq sufficiently removed from a magnetic subband center 
Ec may not be connected to one another through equipo- 
tential islands. Thus if \Eo — Ec\ > W, our quasi-classical 
argument that neglects tunneling among distinct equipo- 
tentials implies a divergent localization length for finite 
edge separation S. The inclusion of quantum tunneling 
would restore the finite value of the localization length ex- 
pected at finite S, but the spatial attenuation of quantum 
scattering between edge states as a function of their sep- 
aration S is then governed by the magnetic length. The 
resulting localization lengths are of a completely different 
order of magnitude than those we find for edge states con- 
nected by the disorder-generated quasi-classical equipo- 
tentials. Indeed, for sufficiently small W, as the Fermi en- 
ergy Eq crosses into the gap \Eo — Ec\ > W the computed 
localization length increases abruptly by many orders of 
magnitude, to values beyond our numerical control. 

Lee and Chalker have proposed a variant of the orig- 
inal one-channel network model to study localization in a 
random-flux background. In order to maintain a vanish- 
ing a-xy, they stipulated that each node describe equiva- 



lent tunneling for the two channels of opposite chirality, 
so that the node is represented by scattering parameters 
6*1, 6*2 satisfying sinh(6'i) = 1/ sinh(02)- They obtain no 
extended state in this model. 

Within the framework of edge-edge scattering de- 
scribed in this paper, their conclusion is readily under- 
stood: no configuration of their nodes can be equivalent 
to a pair of edge states on opposite sides of a long flux 
barrier. Rather, backscattering can occur at every node 
with amplitude equal to the forward scattering that oc- 
curs at that node. In contrast, our backscattering am- 
plitude is exponentially less than our forward scattering 
amplitude, for the entire length of the flux barrier. In a 
random background, this length is set by the quenched 
flux disorder, i.e. the classical path. 

This overall picture supports the general notion, men- 
tioned in our introduction, that the classical paths 
(equipotentials) make up the dominant contribution to 
transport, and that tunneling among equipotentials may 
be viewed as a perturbation of the classical background 

This picture is also consistent with unpublished nu- 
merical studies of one-dimensional defects in Chalker lat- 
tices |l^. In Chalker networks, local isotropy is main- 
tained by an antiferromagnetic choice of node parameter 
6 on alternate sublattices. One can create edges in a 
network by introducing defects in the antiferromagnetic 
ordering, forming domain boundaries across which the 
sublattice alternation is displaced by a half unit-cell. It 
is found that pairs of linear defects (comparable to the 
edges we study here) in a one- channel Chalker model dis- 
play an exponential increase of localization length with 
increasing separation. Two-channel Chalker models dis- 
play a richer behavior. Nodes of the form described 
above (sinh(0i) = l/sinh(6'2)) never exhibit an extended 
regime, whereas nodes satisfying 9 = 9i = 62 with 
sinh(6') = s,l/s on the alternate sublattices |12| show 
decreasing localization lengths as a function of system 
width M for s close enough to 1, and exponentially 
increasing localization lengths as a function of M for 
s > Sc > 1. Details will be given elsewhere ||ll| ]. 

In summary, we have described numerical studies of 
tunneling between one-dimensional edge states created 
by boundaries between regions of differing potential 
and/or flux. These systems display distinct regimes as 
a function of energy and disorder; we have character- 
ized the extended regime, in which localization lengths 
increase exponentially with the spatial separation of the 
edges, and we have offered a quasi-classical explanation 
of this behavior. The fundamental contribution of the 
classical paths, whether formed by magnetic-flux or po- 
tential edges, suggests to us that extended states ought 
to exist in a random-flux background. 
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